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SELECTION  OP  A  DELAY  LINE  MODEL 


by 

rtufu3  Oldenburger^ 
and 

Lawrence  L.  Hoberock2 
Purdue  University 


A  mathematical  model  of  &  linear  system  can  be  derived  using  an 
approximation  of  the  convolution  integral.  This  paper  selects  the 
model  such  that  responses  of  the  model  to  commonly  occurring  inputs 
are  closest  to  corresponding  responses  of  the  system.  The  transfer 
function  of  the  model  ie  the  product  of  the  system  transfer  function  by 
a  linear  combination  of  two  delay  terms  divided  by  an  infinite  product. 
If  the  model  delay  time  is  small  the  infinite  product  may  be  replaced 
by  1,  and  thue  may  be  dropped.  The  delay  time  and  the  number  of  delay 
elements  are  selected  such  that  the  responses  of  the  simplified 
model  are  closest  to  corresponding  responses  of  the  system.  The 
validity  of  the  simplification  is  investigated  for  various  inputs 
by  comparing  the  responses  of  the  simplified  model  with  those  of  the 
exact  model.  It  is  found  for  several  types  of  commonly  occurring 
inputs  that  the  number  of  delay  elements  should  be  chosen  as  large 
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as  physically  possible.  The  result#  show  that  the  value  of  the 
delay  time  should  be  selected  as  a  Ainction  of  the  number  of  delay 
elements  and  of  the  system  bandwidth.  It  is  further  shown  that 
this  function  is  the  same  for  each  of  the  inputs. 
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INTRODUCTION 


To  design  a  controller  for  a  physical  plant,  the  system  must  be 
identified  by  some  method.  An  analytical  representation  of  the  system 
which  may  be  used  to  find  the  output  of  the  system  from  its  input  is 
often  desired,  Truxal*  presents  techniques  for  evaluating  the  differ¬ 
ential  equation,  the  frequency  response  and  the  impulse  response  of 
linear,  time-invariant  systems. 

2 

In  1957  Goodman  and  Res wick  presented  an  application  of  the  theory 
3  A 

developed  by  Tustin  ,  Lewis  and  others  as  a  means  of  identification  for 
linear,  time-invariant  systems.  A  mathematical  model  of  the  system  is 
derived  by  approximating  the  convolution  integral  with  a  weighted  time 
•  series  having  a  delay  interval  T.  The  smaller  the  value  of  T  and  the 

larger  the  number  of  elements  retained  in  the  time  series,  the  better  is 
the  approximation.  Given  a  system  with  unknown  characteristics  the  im¬ 
pulse  response  may  be  approximated  from  the  time  series  as  closely  as 
desired. 

2 

The  delay  line  synthesizer  of  Goodman  and  Keswick  (See  Figure  1.) 
constitutes  a  physical  model  which  uses  only  a  finite  number  of  elements 
in  the  time  series.  The  weighting  factors  are  determined  by  adjusting 

the  synthesizer  settings  until  an  optimum  fit  is  obtained  between  the 
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model  and  system  responses  to  normal  operating  inputs,  Chang  continued 
the  work  in  applying  the  delay  line  synthesizer. 


2 


In  this  paper  the  weighted  time  series  is  used  to  derive  a  mathe¬ 
matical  model  by  keeping  only  a  physically  realizable  number  of  elements 
in  the  series.  The  key  to  this  study  is  the  use  of  an  infinite  product- 
employed  by  R,  Oldenburger  and  R.  E.  Goodson^  for  transcendental  func¬ 
tions  arising  in  the  analysis  of  fluid  flow  through  pipes.  The  Laplace 
transform  of  the  terms  in  the  series  is  taken  and  a  transfer  function 
obtained  by  writing  the  transformation  in  closed  form.  By  choosing  a 
small  value  of  delay  time  the  infinite  product  in  the  transfer  function 
may  be  neglected  and  a  simplified  model  representation  obtained.  Hiis 
approximate  model  representation  is  used  to  compare  system  and  model 
responses  to  various  inputs.  From  this  comparison  the  number  of  delay 
elements  in  the  model  and  the  value  of  the  delay  time  are  selected  such 
that  the  responses  of  the  model  are  closest  to  corresponding  responses 
of  the  system.  The  exact  model  responses  are  compared  with  the  approxi¬ 
mate  responses  to  show  that  the  simplification  is  justified. 

In  what  follows  t  denotes  time. 
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1. 


DERIVATION  OF  THE  DELAY  LINE  MODEL 


Let  g(t)  represent  the  impulse  response,  or  weighting  function,  of 

a  linear  system.  The  response  c(t)  of  the  system  to  an  input  r(t)  may 
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be  determined  by  the  convolution  integral  aB  in 


c(t) 


g(  <r  )  r(t  -  a)  d  <r  ,  t  ,§  0 


(1) 


where 


r(t  -  or)  -  0  ,  t  <  «r  .  (2) 


Let  c^(t)  be  defined  by 

Cj(t)  «  [g(T)  r(t  -  T)  ♦  g(2T)  r(t  -  2T)  *  .  .  . 
♦  g(n  T)  r(t  -  n  T)]t 


(3) 


where  T  is  an  increment  of  the  variable  e*  and  n  is  a  positive  integer. 
If  T  ia  small  and  n  is  large  enough  we  may  write 


c(t)  «*  c^t)  (4) 

where  means  "is  approximately  equal  to."  We  designate  c-^(t)  as  the 
response  of  the  delay  line  model  (DIM)  of  the  system.  The  Laplace 
transform  of  a  function  of  time  f(t),  that  is  jf(t)]  >  is  denotod 

by  F(s),  vherc  s  is  the  Laplace  variable.  From  Equation  (l)  it  may  be 
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shown  that 
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C(e)  R(s)  G(s)  .  (5) 

Sine*  r  is  defined  by  Equation  (2)  we  have® 

f[r(t-nT)]  -  R(s)  e”nTs  .  (6) 

Thus,  taking  the  Laplace  transform  of  the  terms  in  Equation  (3)  yields 
CjU)  -  R(s)  G^s)  (7) 

where 

^(s)  -  [g(T)  e“T“  ♦  g(2T)  e"2T®  ♦  .  .  .  ♦  g(n  T)  e“nTa]T  .  (8) 

I^t  G^(s)  be  defined  as  the  DIM  transfer  function.  Suppose  that  a  linear 
system  has  an  unknown  transfer  function  G(s),  If  the  impulse  response 
g(t)  is  measured  at  times  T,  2T,  .  .  nT,  then  G^Cs )  is  determined. 

The  smaller  the  delay  time  T  and  the  larger  the  number  of  delay  elements 
n,  the  better  is  the  approximation  in  Relation  (3),  or  the  better  does 
G^(s)  approximate  G(s). 

There  must  exist  a  lower  limit  for  T  and  an  upper  bound  for  n  if 
the  DIM  is  to  be  physically  realisable.  Knowledge  of  how  the  difference 
between  system  and  model  responses  depends  on  n  and  T  is  desired.  These 
parameters  are  to  be  selected  such  that  the  model  fulfills  nominal 
engineering  requirements. 
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FIRST  ORDER  SYSTEM 


Consider  the  system  described  by  the  transfer  function 


Employing  the  inverse  Laplace  transform  we  obtain  the  wei#>ting  function 


Hence  G^(s)  in  Equation  (8)  becomes 

01(e)  -  a  T  [e-T<a+8>  ♦  e“2T<a+e>  ♦  .  .  .  ♦  e-nT<a+s)]  .  (ll) 


Using  the  transformation 


z  »  T(a  +  s) 


we  write  Equation  (ll)  in  the  closed  form 


1  T  S  eiz  _  e“^z 


It  is  desired  to  expand  the  denominator  of  the  term  on  the  right  of 
equation  (13)  as  an  infinite  product.  Let  a  function  Z(z)  have  simple 
zeros  at  the  points  z^,  zg*  *3>  •  •  •  >  where 


Lim 

k  •+  oo 
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these  being  the  only  zeros  of  Z(z)  and  i 0.  If  Z(y.)  is  analytic 

9 

for  all  values  of  z,  it  may  be  shown  that. 


where  exp  £  j  denotes  raising  the  base  "e"  to  the  power  indicated  in 
the  brackets.  The  zeros  of  the  denominator  above  are  given  by 

^  -  2klTj  ,  k  -  0,  i  1,  i  2,  .  .  .  (15) 

where  j  *  \/-l.  Hence  the  conditione  for  expansion  of  the  denominator 

of  Equation  (13)  by  Equation  (14)  are  not  satisfied  due  to  the  root 
as  *  0.  This  problem  is  solved  by  writing 

>  -  -  *  ~Ve"*g  .  <16> 

The  bracketed  expression  now  satisfies  the  conditions  for  expansion,  the 
zeros,  zr,  being  given  by 

zr  «  2rTTj,  r  ■  ♦  1,  +  2,  .  .  .  _  (17) 

Thus  by  Equations  (14)  and  (16) 

.  ,-i*  .  z  [l  *  ^*2-]  [l  -  ~pr]  ...  .  (18) 

Substituting  the  expression  for  z  in  Equation  (13)  yields 


V*1  ‘[rh]  [I] 


(19) 
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inhere 

1 _ 

r  .  T2(a^:  | 

l  a2  *2  J 

H(s)  -  [e-^aT  e-£Ts  -  e-a(n^)T  e~<n+i)TsJ  (23.) 

and  the  index  m  is  a  non-negstive  integer.  Figure  2  presents  a  block 

diagram  of  G^(s)  in  Equation  (19). 

* 

Let  a  new  function  C^(s)  be  given  by 


ci(s)  •  r(s)[i4t1  fl" 


(22) 


Then  from  Equations  (7)  and  (19)  it  follows  that 
C^Cs)  -  cj(s)  H(s) 


(23) 


Substituting  the  expression  for  H(u)  and  employing  the  inverse  Laplace 
transform  gives 

cx(t)  -  e^aTCi[t-iT]  u  [t-jT] 


-  e“a<n+*)T  cj  [t  -  (n  ♦  J)  T  ]  u  [t  -  (n  ♦  *)  t] 


(24) 


where  c^(t)  is  the  inverse  Laplace  of  C^(s)  and  the  function  u  is  defined 
by 

'  0,  t  <  ti 


u(t  -  t]_)  *  < 


1.  t  >  tl 


tl  >  0 


.  (25) 
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3.  APPROXIMATION  OF  THE  DIM  TRANSFER  FUNCTION 

Let  T  be  small  enough  such  that  the  product  aT  is  much  less  than 
one  (aT  «  1).  Then  a  reasonable  approximation  of  G^(s)  results  when 
the  infinite  product  X  is  neglected.  From  Equation  (22) 

c£(a)  ft*  C(s)  ,  aT  «  1  .  (26) 

We  introduce  C^(s),  where 

C1A(s)  -  C(e)  H(s)  .  (27) 

Hence  from  Equation  (23) 

C^(e)  &  Cj^(e)  ,  aT  «  1  .  (28) 

I^t  the  approximate  model  response  be  denoted  by  c^A(t),  which  is  the 
inverse  Laplace  transform  of  Cia(s),  Substituting  the  expression  for 
H(s)  in  Equation  (21)  and  employing  the  inverse  Laplace  transform  we 
obtain 

c1A(t)  -  e^aT  e[t  -  £t]  u[t  -  Jt] 

(29) 

-  e“a(n**)T  c[t  -  (n  ♦  £)t]  u  [t  -  (n  *  £)t]  . 

To  differentiate  from  the  approximate  model  response  we  henceforth  refer 
to  c^(t)  as  the  exact  model  response.  The  problem  of  selecting  n  and  T 
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is  simplified  using  Relation  (28).  In  the  following  study  we  select  the 
nodel  parameters  and  show  that  employing  this  relation  is  reasonable  for 
five  conmonly  occurring  inputs  to  the  system. 
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STEP  INPUT 


Consider  a  step  input  to  the  model  snd  the  system  given  by 

r(t)  -  u(t)  .  (30) 

Using  the  Laplace  transform,  from  Equation  ($)  we  obtain  the  system 
response  as 

c(t)  -  [l  -  e"4t,j  u(t)  .  (31) 

From  Equation  (29)  the  approximate  model  response  becomes 


o 

1 

's 

o 

,  0  i  t  <  & 

(32a) 

cu(t)  -  e-*aT  -  e-at 

,  JT  ^  t  <  (n  4  £)T 

(32b) 

cu(t)  -  e-**T  [l  -  e-*nT] 

i  t  £  (n  4  £)T  , 

(32c) 

Selection  of  Model  Parameters 

By  Relations  (32)  the  approximate  model  response  varies  with  time 
only  for  £  t  <  (n  4  ^)T.  For  t  outside  these  limits  the  response 
of  the  model  fails  to  follow  that  of  the  system.  Figures  3  through  5 
show  the  system  and  model  responses  to  a  step  input  of  magnitude  K  with 
unity  bandwidth  (a  *  l)  and  various  values  of  T  and  n.  The  responses  are 
normalized  by  dividing  each  term  by  K.  The  model  response  is  more  satis¬ 
factory  for  large  n  and  small  T,  as  for  example,  n  *  300  and  T  *  0.01  in 


SYSTEM  RESPONSE 


APPROX.  MODEL 
RESPONSE  - 


System  and  Approximate  Model  Responses 
Step  Input  for  T  “  0.1 
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Figure  3.  However,  a  delay  line  model  with  such  parameters  would  be 
unpractical.  Goodman  and  Reswick  use  n  *  20  in  the  delay  line  synthe- 
sizer  .  We  see  from  the  response  plots  mentioned  above  that  for  any  n 
the  final  value  of  ciA(t)  has  a  maximum  which  is  dependent  on  T.  Ws 
designate  T0  as  that  value  of  T  for  which  the  maximum  final  value  occurs. 
As  T  -*■  Te  from  values  less  than  T0,  the  final  value  of  c^A(t)  increases, 
but  the  model  follows  the  system  lees  closely  for  the  range 
£T  £  t  <  (n  ♦  £)T.  For  finite  non-zero  values  of  n  and  T  the  model 
output  is  always  less  than  that  of  the  system.  We  choose  to  select  n 
and  T  to  maximize  the  final  value  of  C]_^(t).  Since  for  any  T  the  best  n 
to  accomplish  this  is  infinitely  large,  we  first  select  n  as  large  as 
reasonably  possible  for  applications.  The  time  T  is  then  chosen  to 
maximize  the  final  value,  7,  of  the  approximate  model  response,  where 
from  Equation  (32c) 

V  -  e^aT  [l  -  •_anT]  •  (33) 

Let  ne  be  the  value  of  n  selected  for  the  model.  The  maximum  value  of  V 
occurs  when  T  -  T0,  where 


Ir(2np  +  1) 

a  no 


(34) 


The  maximum  final  value,  Vm,  of  0^(1 )  i*  obtained  by  replacing  T  in 
Equation  (33)  by  the  expression  for  T0,  giving 

Vjj  -  (2no  *  -  (2no  ♦  if1  ]  .  (35) 

If  n©  ■  20,  the  value  used  by  Goodman  and  Reswick,  and  a  -  1,  for 
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simplicity,  Equations  (34)  and  (35)  yield 
T0  -  0.1857  ,  Vm  -  0.8891 

Figure  6  indicates  system  and  approximate  model  normalized  responses  for 
a  step  input  of  magnitude  K,  with  a  ■  1  and  T  »  T0  at  no  ■  20.  It  is 
easily  verified  from  Equations  (28)  and  (29)  that 


lira 

n0-*<» 


To 


0 


(36) 


Lim 
n0-»  oo 


vm 


i 


(37) 


Thus  better  model  response,  using  values  of  T0,  is  obtained  by  increasing 
n. 


Comparison  of  the  Approximate  Model  with  the  Exact  Model 
We  have  selected  the  model  parameters  using  an  approximation  which 
neglects  the  infinite  product  T8  in  the  expression  for  Gi(s).  To 
determine  the  effects  of  using  this  approximation  we  replace  R(s)  in 
Equation  (22)  by  the  Laplace  transform  of  the  step  input  to  obtain 


*/  v 

Cj^s) 


a 

s  ♦  a 


(38) 


Expanding  this  expression  for  ( s )  in  partial  fractions  and  employing 
the  inverse  Laplace  transform,  we  have 


c*(t) 


1L 


-at 


+  Yl 

m=l 


sin 


2TTm 


t  -  0, 


m 


(39) 


20 


ci(t) 


1  -  e" 


anT  ,'iaT 


t  £  (n  +  £)T  .  (46c) 


Comparison  of  Equations  (46)  with  Equations  (32)  indicates  that  both 
c^(t)  and  c^A(t)  are  zero  for  0  £  t  <  j(T.  For  the  second  range, 

Kt  <  (n  ♦  J)T,  the  effect  of  the  approximation  is  to  omit  the 
attenuating  infinite  product  1T"t  in  the  first  term  and  the  infinite 
sinusoidal  series  in  the  second  term  of  the  exact  model  response.  For 
t  j>  (n  *  £)T  the  approximation  omits  Tt  in  the  exact  model  final 
value.  We  note  that  Hr  -+  1  as  T  -»  0.  Further,  Tn  -  (-iK1  2. 
For  small  a  T  we  have  «m  *  **  °*  Taking  ■  0, 

dm  "  the  summation  in  Equation  (39)  i»T3  between  zero  and 

•&T.  depending  on  the  value  of  t.  The  same  is  true  for  the  summation  in 
Equation  (46b).  These  results  and  numerical  studies  indicate  that  if 
aT  <  £  °ne  can  neglect  the  summation  and  hence  the  product  TTS. 

To  evaluate  graphically  the  effects  of  the  approximation  we  employ 
Equation  (3)  for  the  exact  model  response.  Let  the  inverse  Laplace 
transform  f(t)  of  a  function  F(e)  be  denoted  by  J[_1  [?(»)] .  Since8 


r  [<-f>  (•■nT*)]  *«<»-»*> . 


Equations  (3)  and  (10)  are  used  to  give 


a  T  [e~aT  u(t  -  T)  «■  e_2aT  u(t  -  2T)  +  .  .  . 
e“naT  u(t  -  nT)] 


Figure  7  shows  the  system  and  exact  model  normalized  responses  for  a 
step  input  of  magnitude  K,  with  a  ■  1,  T  -  T0  and  *  20.  Comparing 

these  plots  with  those  in  Figures  3  through  6  indicates  that  the  most  sig¬ 
nificant  effect  of  the  approximation  is  the  smoothing  of  the  "staircase" 


System  and  Exact  Model  Responses  to  a  Step 
Input  for  T  ■  T0  at  n0  ■  20 
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response  of  the  exact  model.  A  farther  effect  of  the  approximation  is 
to  decrease  the  time  at  which  the  time-variant  model  output  begins  from 
T  to  £T.  The  time  at  which  it  ends  is  increased  from  n  T  to  (n  ♦  £)T, 

This  as  well  as  the  smoothing  is  attributed  to  omitting  the  sinusoidal 
series  in  the  approximate  model  response. 

We  wish  to  determine  whether  the  results  in  selecting  the  approxi¬ 
mate  model  are  applicable  to  the  exact  model.  It  is  recalled  that  n  was 
dioeen  as  large  as  reasonably  possible  for  applications,  lhe  delay  time 
T  was  then  selected  such  that  the  final  value  of  the  approximate  model 
response  given  by  Equation  (33)  was  maximized.  From  Equation  (46c)  the 
exact  model  final  value,. V^,  is  given  by 

71  -  [Hr]  [l  -  •“anTj  •■*aT  .  (49) 

An  attempt  to  maximize  thia  expression  for  by  differentiating  with 

respect  to  T  would  be  fruitless.  We  choose  to  select  n  and  T  for  the 

exact  model  the  same  as  for  the  approximate  model  provided  that  the 

value  of  IT  is  approximately  unity.  ThiB  means  that  the  product  a  T 

must  be  small.  It  is  desired  to  know  how  small.  In  his  discussion  of 

2 

the  paper  by  Goodman  and  Reswick  ,  C.  M.  Chang  notes  that  if  the  delay 
time  T  is  greater  than  ■— ,  where  a  is  the  system  bandwidth,  the  fre¬ 
quency  response,  particularly  in  the  high  frequency  region,  obtained  on 

the  delay  line  synthesizer  is  questionable.  This  is  derived  from  the 

10 

fundamental  theorem  of  sampling  which  states  that  a  signal  is  completely 
determined  by  values  of  the  signal  (samples)  taken  at  a  series  of  in¬ 
stants  separated  by  T^  -  £a,  where  T^  is  the  sampling  period.  From  these 
considerations,  for  the  remainder  of  this  work  it  will  be  assumed  that 
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0  <  a  T  <  i 


(50) 


Lst  new  quantities  and  be  defined  by 

^  -  nr^ 


(51) 


\  -  £  tn  [l  ♦  \  ]  .  (52) 

m*l 

Then  for  the  denominator  of  the  expression  for  Tl^  in  Equation  (4.0)  we 

have 


¥  [l  ♦  *m]  " 


But  it  is  seen  that 


k.  <  E  W  .  i. » o 


Renee  from  the  relation 

y  i  .E 

•  *  m  2  Z 

m-1  a  6 


we  may  write 


Vfc  <  ~  ,  0<a2l2<i 


(53) 


(54) 


(55) 


(56) 


This  gives 


“5E"  >  0.9900 


(57) 


and  we  conclude  from  Equation  (40 )  that 
0.9900  <1^.  <  1 

Consequently  It  Is  safe  to  assume  that  the  results  in  selecting  the 
approximate  model  are  reasonably  applicable  to  the  exact  model. 
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5. 


RAMP  INPUT 


The  ramp  signal  given  by 

r(t)  -  t  u(t)  (59) 

ie  the  second  input  to  be  considered.  Using  transform  techniques  we 
obtain  from  Equation  (5) 

e(t)  -  [t-i.  ♦  — #“atl  u(t)  .  (60) 

•■a  a  J 

To  derive  the  approximate  model  response  we  substitute  this  expression 
for  c(t)  in  Equation  (29),  giving 

cu(t)  -  0  ,  0  £  t  <  US  (61a) 

cu(t)  -  e4aT  [t  -  Jt  -  -i-1  ♦  i  e“At  , 

(61b) 

it  £  t  <  (n  *  i)T 

cu(t)  -  [e^aT-e-a<n^)T]t 

-  e-i»T  hr  *  J  -  e-mT  [(n  ♦  J)T  ♦  i]  }, 

1  J  (61c) 


t  g  (n  ♦  J)T 
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Selection  of  Modal  Parameters 

Comparison  of  the  approximate  model  response  with  the  system  res¬ 
ponse  indicates  that  for  small  T  the  response  of  the  model  follows  that 
of  the  system  reasonably  well  for  the  second  range,  iT  <  t  <  (n  *  J)T. 
For  the  first  period,  0  £  t  <  gT,  the  approximate  model  output  is 
zero,  as  in  the  step  input  case,  while  for  the  last  range, 
t  >  (n  *  £)T,  a  lag  in  slope  as  well  as  position  appears  in  the  response 
of  the  model.  Figures  8  and  9  show  the  system  and  approximate  model 
normalized  responses  for  a  ramp  input  of  magnitude  K  with  various  values 
of  T  and  n  and  unity  bandwidth.  Again,  large  n  and  small  T  give  more 
satisfactory  model  response,  however  these  values  are  limited  in  practi¬ 
cal  applications.  For  email  T  the  approximate  model  response  follows 
that  of  the  system  satisfactorily  during  the  second  period.  For 
t  £  (r.  ♦  £)T  the  approximate  model  response  exhibits  a  constant  final 
slope,  S,  where  from  Equation  (6ie) 

S  .  e“*aT  [l  -  e~anT]  .  (62) 

The  system  response,  however,  has  a  constantly  increasing  slope  approach¬ 
ing  unity  as  t-*  oo  ,  The  expression  for  S  in  Equation  (62)  is  identical 
with  that  for  V  in  Equation  (33).  Hence  S  has  a  maximum  value  occurring 
at  T  ■  T0.  As  T  •+  Te  from  values  less  than  T0,  the  slope  S  increases; 
however,  the  response  c^(t)  follows  the  system  response  c(t)  less  close¬ 
ly  for  £  t  <  (n  *  £)T.  Again,  selecting  the  model  parameters  pre¬ 
sents  a  conflict  between  maximizing  S  and  improving  the  model  response 
for  the  second  range  of  t.  We  dioose  to  select  n  and  T  to  maximize  the 
final  slope  of  the  sp  proximate  model  response.  Thus  n  and  T  are  selected 
the  same  as  in  the  step  input  case. 


System  and  Approximate  Model  Responses  to 
Input  for  T  =  T0  at  Ho  “  20 
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Comparison  of  the  Approximate  Model  with  the  Exact  Model 
It  is  desired  to  determine  the  effects  of  the  approximation  in 
which  "fljj  in  the  expression  for  °1<  s)  in  Equation  (19)  was  omitted. 
Replacing  R(s)  in  Equation  (22)  by  the  Laplace  transform  of  the  ramp  in¬ 
put  yields 


cj(s) 


[■?]  [rh-]  R 


(63) 


By  expandir^  this  expression  for  cj(s)  in  partial  fractions  and  employ¬ 
ing  the  inverse  Laplace  transform  we  obtain 


(64) 


where  I,  ,  Ajj,  and  6^  are  given  in  Equations  (40),  (41),  and  (42) 
respectively.  From  Equations  (24)  and  (64)  it  follows  that 


Ci(t)  ■  0  , 


0  5  t  <  jfl 


(65a) 


ei(t)  -  [  1Tt][.^T  (t  -  4T)]  -  ±  «^aT 


♦  e 


-at 


(65b) 


£T  S  t  <  (n  ♦  £)T 
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(65c) 


Comparison  of  Equations  (65)  with  Equations  (61 )  indicates  again  that  the 
effect  of  the  approximation  is  the  absence  of  the  attenuating  product 
Tt  from  the  sp  proximate  model  response.  Also  missing  is  the  infinite 
sine  series  which  causes  the  discontinuities  in  the  exact  model  response. 

To  illustrate  graphically  the  effects  of  the  approximation,  the 
expressions  for  r(t)  and  g(t)  in  Equations  (59)  and  (10)  are  used  in 
Equation  (3),  giving 


ci(t)  *  a  t  [e“a^  (t  -  T)  u(t  -  T)  ♦  e-2aT  (t  -  2T)  u(t 
♦  .  .  .  +  (t  -  n  T)  u(t  -  n  T)]  . 


2T) 


(66) 


Figures  10  and  11  show  the  system  and  exact  model  normalized  responses 
for  a  ramp  input  of  magnitude  K  with  a  ■  1  and  several  values  of  n  and 
T,  It  is  recalled  that  the  approximation  for  the  step  input  resulted  in 
the  smoothing  of  the  "staircase"  exact  model  response.  For  the  ramp  in¬ 
put,  however,  comparison  of,  for  example,  the  responses  in  Figure  8  with 
those  in  Figure  10,  shows  a  smoothing  of  the  "staircase  slope"  of  the 
exact  model  response.  A  further  effect  of  the  approximation  is  to  de¬ 
crease  the  time  at  which  the  model  response  begins  from  T  to  ^T.  Also, 
the  time  at  which  the  constant-slope  portion  of  the  model  response 


System  and  Exact  Model  Responses  to  a  Ramp 
Input  for  T  -  0.1 


32 


System  and  Exact  Model  Responses  to 
Ramp  Input  for  T  *=  1.0 
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begins  is  increased  from  nj  to  (n  +  £)T.  This  and  the  smoothing  are 
again  attributed  to  the  absence  of  the  infinite  sine  series  in  the  ap¬ 
proximate  model  response. 

In  Figures  10  and  11  the  staircase  nature  of  the  slope  for  the  exact 
model  response  is  indicated  by  dashes.  These  dashes  are  straight  line 
segments  and  are  connected.  The  connections  are  left  out  to  make  the 
changes  in  slope  clear  to  the  reader.  Thus  in  Figure  10  each  dash  lasts 
0.1  unit  of  time,  tee  dash  runs  from  t  ■  0.1  to  t  “  0.2,  a  second  from 
t  *  0.2  to  t  =  0.3,  etc.  In  the  case  of  the  n  =  20  curve,  which  starts 
at  t  *■  0.1,  the  slope  changes  every  0.1  unit  of  time  until  t  «•  2.0.  After 
t  *  2.0  the  n  «  20  curve  has  a  constant  slope. 

As  in  the  step  input  study  the  concern  is  whether  the  results  in 
selecting  the  approximate  model  are  applicable  to  the  exact  model.  We 
chose  to  first  select  n  as  large  as  possible  for  practical  applications. 
The  time  T  was  then  selected  to  maximize  the  final  slope  S  of  the  ap¬ 
proximate  model  response.  Biis  gave  the  same  result  for  To  as  in  the 
step  input  study.  From  Equations  (65c)  and  (62)  the  final  slope,  Si, 
of  the  exact  model  response  is  given  by 


Hence  by  Relation  (58)  we  conclude  that  maximizing  S  rather  than  Si  is 
a  reasonable  method  for  selecting  the  value  of  the  delay  time. 
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IMPULSE  INPUT 


The  third  input  to  be  considered  is  given  by 


r(t)  -  $(t) 


where  the  unit  impulse,  6  (t),  is  defined  by 


6  (t  -  tj^)  -  0  , 


rt+t. 

Js( t  -  tj.) 


d  t  *  1  , 


t  h  *3 


6  >  0  J  . 


Since  the  input  is  an  impulse,  the  weighting  function  in  Equation  (10) 
is  used  to  give  the  system  response 


c(t)  ■  a  e"at  u(t) 


Hence  by  Equation  (29) 


Cjj^(t)  *0  , 


oj|t  <  £t 


ou(t)  -  a  e’at  ,  &  f  t  <  (n  ♦  £)T 


Ci«(t)  ■  0  , 


t  >  (n  *  in 


(71a) 


(71c) 


35 


Selection  of  Model  Parameters 

For  the  period  J  ^  t  <  (n  +  |)T,  the  approximate  model  response 
ia  identical  with  the  response  of  the  system.  Figure  12  indict tee  the 
system  and  approximate  model  normalized  responses  for  an  impulse  of 
strength  K  with  a  ■  1  and  various  values  of  n.  Only  one  value  of  T  is 
used  due  to  the  simplicity  of  the  response  plots.  Good  model  response 

m 

results  when  the  lower  limit  of  time  t  in  Equation  (71b)  is  small  and 
the  upper  limit  (n  +  £)T  is  large.  That  is,  the  smaller  the  value  of  T 
and  the  larger  the  value  of  n,  the  better  is  the  model  response.  Prac¬ 
tical  values  of  these  parameters  cause  conflicts  in  model  selection.  For 
any  T  the  best  n  is  as  large  as  possible  in  physical  applications.  On 
the  other  hand,  the  problem  of  selecting  the  best  T  for  any  n  is  that  if 
T  is  chosen  too  small  the  interval 


(n  +  £)T  -  &T  -  n  T 


of  good  model  response  is  too  small.  If  T  is  selected  too  large,  the 
starting  time  t8  *  of  the  model  output  is  too  great.  Hence  the  dif¬ 
ficulty  in  selecting  T  for  the  impulse  study  differs  from  that  in  the 
previous  two  cases.  The  criterion  for  selection  is  time  interval  and 
starting  time  of  model  response,  or  graphical  abscissa,  rather  than 
values  of  model  output  or  graphical  ordinate.  We  choose  T  to  maximize 
the  area  under  the  approximate  model  response  curve  as  a  reasonable 
solution  to  the  problem.  The  area  A^  under  this  curve,  found  by  inte¬ 
gration,  is  given  by 


*i  “ 


e-£al  |"i  _  e-anTj 


(72) 


System  and  Approximate  Model  Responses  to  an 
Impulse  Input  for  T  ■  0.1 
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This  expression  is  identical  with  that  for  V  and  S  of  the  previous  two 
cases.  Hence  the  value  Te  of  the  delay  time  T  to  give  maximum  A^  is 
given  by  Equation  (34).  The  expression  for  maximum  area  is  the  same  as 
that  for  Vm  in  Equation  (35).  from  Equations  (36)  and  (37),  increasing 
flowers  T0  and  increases  the  maximum  area  as  before. 

Comparison  of  the  Approximate  Model  with  the  Exact  Model 
To  determine  wh ether  using  the  approximation  is  a  reasonable  method 
of  selecting  n  and  T,  the  exact  model  response  must  be  determined.  Since 
R(s)  *  1  for  the  impulse  input,  Equation  (22)  becomes 

’  [r4r][T.]  •  to) 

Using  partial  fraction  expansion  and  the  inverse  Laplace  transform  we 
obtain  from  Equation  (73) 


c,(t) 


.-at 


-  «* 

{**  ^ 

l  m-1 


Am  sin 


iJjl.  t  ♦  e 

iy  m 


} 


(74) 


Thus  by  Equation  (24) 


L(t)  -  0  , 


cx(t) 


0  £  t  <  JT  (75a) 


[.  •  V.  K.  'In  [4^  t  •  -  „T 


e"at  ia  *  Z_,  sin 

m-1 


[  T 


(75b) 


jiT  <  t  <  (n  +  £)T 


=l(t) 


t  2  (n  ♦  £)T 


0 


> 


(75c) 
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Unlike  the  responses  in  the  step  and  ramp  input  studies,  the  infinite 
product  IT  does  net  appear  in  the  expression  for  c^(t)  for  the  im¬ 
pulse  input.  Hence  the  effect  of  omitting  TTa  from  G^s)  in  the 
simplification  technique  is  the  absence  of  only  the  infinite  sine 
series  from  Equations  (75). 

To  facilitate  plotting  e^Ct)  versus  t.  Equation  (3)  is  used  to  give 
for  the  impulse  input 


e^(t)  -  a  T  [e"aT  *(t  -  T)  +  *-2»T  S(t  -  2T) 

♦  .  .  .  ♦  e~mT  S  (t  -  n  T>] 


(76) 


The  system  and  exact  model  normalized  responses  for  an  impulse  input  of 
weight  K  are  shown  in  Figure  13,  with  a  °  1,  T  »  0.1,  and  various 
values  of  n.  The  exact  model  response  is  a  series  of  impulses  corres¬ 
ponding  to  a  sampling  of  the  system  response.  Hence  the  effect  of  the 
approximation  ahown  in  Figure  12  ie  that  of  an  ideal  holding  device 
over  a  limited  time  interval.  Omitting  the  Infinite  sine  series  in  the 
approximate  model  response  not  only  changes  the  limits  of  model  response 
from  T  to  £T  and  nT  to  (n  +  |T)T  but  also  removes  the  discontinuities 
of  the  exact  model  response. 

Since  the  ordinates  of  the  response  ci£t)  represent  strengths  of  im¬ 
pulses,  there  is  no  area  under  the  response  curve  to  maximize  in  selecting 
T.  The  problem  of  selection,  however,  is  similar  to  that  in  the  appro¬ 
ximation.  The  interval  nT  -  T  is  desired  large,  which  for  constant 
n  means  increasing  T.  Conversely,  the  starting  time  of  model  response 
and  interval  between  the  impulses  is  desired  small.  The  exact  model 
response  is  a  pure  sampling  of  the  system  response  for  the  time  interval 
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System  and  Exact  Model  Responses  to 
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of  Equation  (75b).  Thus  a  technique  for  selecting  the  delay  time  is  to 
maximize  the  area  under  the  system  response  curve  between  limits  deter¬ 
mined  by  the  exact  model  response,  namely  j^T  <  t  <  (n  ♦  £)T.  But  this 
area  is  exactly  that  given  in  Equation  (?2).  Thus  unlike  the  previous 
two  input  cases,  selection  of  the  delay  time  could  have  been  effected 
without  using  the  approximation. 
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7.  PULSE  INPUT 

The  last  transient  input  to  be  considered  is  the  pulse  shown  in 
Figure  14a.  For  K  ■  1,  where  K  is  the  height  of  the  pulse,  the  input  is 
given  by 

r(t)  -  u(t)  -  u(t  -  T)  (77) 

where  T  is  the  duration  of  the  pulse.  It  is  easily  shown  that  the  sys¬ 
tem  response  to  this  input  is  given  by 

c(t)  -  [l  -  e“at]u(t)-[l  -  e"*^”  u(t  -  r)  .  (78) 

From  Equation  (29)  it  follows  that 


CU(t) 

-  o  , 

0  <  t 

< 

(79a) 

clA^t) 

-  e"*aT  -  e-at  , 

i*  i  t 

<  (T*  jT) 

(79b) 

T  <nT« 

-  X)  m  e-at  # 

.OlA^ 

(T+  £t)  <  t 

<  (n  ♦  i)T 

(79c) 

7*>  n  T  |c!A(t) 

-JaT  -at 

-  e  -  e  , 

i*  <  t 

<  (n  ♦  J)T 

(79d) 

°lA(t) 

-  e^aT  [l  -  e-anT]  - 

[e-i®T  -  e-a^" 

r^]  u(t  -  t 

-4») 

(n  ♦  J)T  <  t  <  7  ♦  (n  +  £)T 

(79e) 


ciA(t)  -  0 
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t  fc  T*  (n  ♦  J)T  .  (79f) 

Selection  of  Model  Parameters 
Comparing  Equations  (79)  with  Equation  (78)  reveals  again  that  the 
smaller  the  value  of  T  and  the  larger  the  value  of  n,  the  better  is  the 
approximate  response.  Figure  14b  shows  system  and  approximate  model 
normalized  response  for  a  pulse  input  of  magnitude  K  with  a  *  1  and 
values  of  T,  n,  and  't  as  shown.  The  output  of  the  model  is  equal  to 
that  of  the  system  for  (£r  ♦T')  it  £  (n  ♦  £)T.  For  t  outside  these 
limits  the  value  of  c^A(t)  is  always  less  than  c^(t).  For  any  T  the 
best  n  is  as  large  as  practical  in  physical  applications.  If,  on  the 
other  hand,  n  is  held  constant  and  T  is  reduced  to  improve  the  response 
of  the  approximate  model  for  it  <  X  *  JT,  the  response  for 
t  £  (n  ♦  i)T  is  less  desirable.  If  T  is  increased  the  reverse  is  true. 
As  a  reasonable  compromise  to  the  conflict  in  selecting  the  best  delay 
time,  pick  n  as  large  as  possible.  Then  choose  T  to  maximize  the  area 
under  the  approximate  model  response  curve.  By  integrating  it  may  be 
shown  from  Equations  (79)  that  thie  area,  Ap,  is  given  by 

Ap  -  r  [l  -  e~anTj  e-£aT  .  (go) 

With  the  exception  of  the  pulse  duration  ,  the  expression  for  Ap  is 
the  same  as  that  for  ¥,  S,  and  encountered  earlier.  Hence  the  delay 
time  to  give  maximum  area,  Ap^,  is  given  by  Equation  (34)  and 

Apa  -  tVB  (81) 


where  Vm  is  given  in  Equation  (35) 
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Figure  14a  -  Pulse  Input 

Figure  14b  -  System  and  Approximate  Model  Responses  to  Pulse  Input 
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Comparison  o f  the  Approximate  Model  with  the  Exact  Model 
Substituting  the  expression  for  the  impulse  input  and  the  system 
weighting  function  in  Equation  (3)  we  obtain 

c^t)  -at  £e“aT£u(t  -  T)  -  u(t  -  T  -  T)]  u(t  -  T) 

♦  e~2aT  [u(t  -  2T)  -  u(t  -  2T  -  r)]  u(t  -  2T)  (82) 

*  ...  *  e'^  £u(t  -  n  T)  -  u(t  «  n  I  «  '!')]  u(t  -  nT)|. 

By  inspection,  the  area  A]_  under  the  curve  given  by  a  plot  of  cj_(t)  is 
given  by 

A1  '  *  [a  T]  [e-4T  ♦  e"2aT  ♦  •  .  •  ♦  e“naT]  .  (83) 

But  the  product  of  the  bracketed  quantities  in  this  expression  is  easily 
seen  to  be  the  final  value  of  c-^(t)  in  Equation  (48).  From  Equation 
(49)  it  follows  that 

A1  "  t[¥T]  [l  -  e-anT]  e"^aT  .  (84) 

By  Relation  (58)  it  is  concluded  that  choosing  T  to  maximize  Ap  in 
Equation  (80)  rather  than  is  a  reasonable  technique  for  selecting 
the  delay  time. 
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8. 


SINUSOIDAL  INPUT 


Consider  lastly  the  sinusoidal  input  given  by 

r(t)  -  sin  (  (o  t)  u(t)  (85) 


where  to  is  the  signal  frequency.  Let  the  quantities  Q^,  Qg,  and  V 
be  defined  by 


a  q) 

a2  ♦  u>2 

Q2  “  (a2  «■  cj2)jf 

*  *  tan*--1-  (— ) 


(86) 

(87) 

(88) 


Then  for  the  approximation  study  we  follow  procedures  similar  to  those 
in  the  previous  four  cases  to  obtain 

c1A(t)  -  0  ,  0  <  t  c  £T  (89a) 

c,.(t)  ■  Q-.e“at  ♦  Q2e”£aT  sin  [cO  t  -  <f>  -  ^  t]  , 

L  (89b) 


£T  <  t  <  (n  ♦  £)T 
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c^Ct)  -  Q2e”^aT  £sin  (a)t  -  <f  - 

-  e*6nT  sin  (  u)  t  -  *  -  £  u>r  -  a) n  T)]  ,  (89c) 


t  i  (n  ♦  J)I 


Selection  of  Model  Parametera 

The  response  of  the  system  to  the  sinusoidal  input  is  given  by 

c(t)  -  8in  [e<>t  “  t  ]|  u(t)  • 

From  Equation  (89b)  decreasing  the  value  of  T  improves  the  model  response 
during  the  specified  interval  by  increasing  the  value  of  e”^aT  toward 
one  and  decreasing  the  value  of  the  phase  lag  £  to  T  toward  zero.  On  the 
other  hand,  from  Equations  (89b)  and  (89c),  too  small  a  value  for  T  gives 
poor  over-all  model  response.  For  eny  T  the  number  n  should  be  chosen 
as  large  as  possible.  Again  some  compromise  is  needed  in  selecting  T. 

We  compare  c^A(t )  with  c^(t)  for  t  £  (n  ♦  |)T.  Assume  that  the  quantity 
(n  ♦  i)T  is  large  enough  such  that  the  exponential  in  Equation  (90)  may 
be  neglected,  giving 

c(t)  <52  sin  [wt  -  p]  ,  t»-~  •  (91) 

The  first  sine  term  of  Equation  (89c)  contains  a  phase  lag  of  £  u)  T 
radians  from  the  sine  term  of  the  system  response.  The  second  sinusoid 
lags  the  first  by  ui  n  T  radians.  However,  assuming  that 

e-^aT  »  e-a(n^)T  (92) 


the  effect  on  cj_g(t)  of  the  second  sinusoid  of  Equation  (89c)  will  be 
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negligible  regardless  of  phase  shift  a)n  T.  Conversely,  if  the  second 
sine  terra  with  its  attenuating  exponential  is  neglected  completely  the 
parameter  n  will  be  absent  from  the  remaining  terms.  No  method  for 
selecting  T  would  then  be  available.  Hence  as  a  compromise  we  write 


cu(t)  Qjje-^aT  £l  -  e"anT]  sin  £u)t  -  <f  -  , 

t  £  (n  ♦  i  yi  5  <n  ♦  £)T  »  ~ ;  n  T  »  JT 

a 


(93) 


in  which  the  effect  of  the  phase  shift  u>n  T  is  neglected.  Comparison 
of  Relation  (93)  with  Relation  (91)  indicates  that  T  should  be  chosen  to 
minimize  the  shift  T  and  maximize  the  quantity  e“^aT  [l  -  e“an^j. 
However,  since  the  system  is  linear,  the  delay  time  should  not  be  a 
function  of  the  signal  frequency  Hence,  we  select  T  only  to  max~ 

imize  the  attenuating  expression  aoove.  The  input  frequency  must  be 
kept  small  in  order  that  the  model  response  phase  lag  be  small.  Thus 
the  model  parameters  n  and  T  are  selected  as  in  the  previous  four  cases, 
the  delay  time  TQ  being  given  by  Equation  (34).  Figure  15  shows  the 
system  and  approximate  model  normalized  responses  as  given  by  Relations 
(91)  and  (93)  for  a  sinusoidal  input  of  magnitude  K.  The  values  of  the 
parameters  used  are  a  ■  l,u>  ■  1,  n  ■  20,  and  T  »  T0. 


Comparison  of  the  Approximate  Model  with  the  Exact  Model 
To  determine  the  exact  model  response  we  replace  R(s)  in  Equation 
(22)  by  the  Laplace  transform  of  the  sinusoidal  input  to  obtain 


Cj(e) 


u> 

a 

g2  * 

.  s  ♦  a 

(94) 


kS 


Is 


System  and  Approximate  Model  Responses  Given  by 
Relations  91  and  93  to  a  Sinusoidal  Input 
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Let  the  quantity  be  the  absolute  value  of  B^  where 


**  "  *  *  3“  I  j 


r  T2(a  ♦ 

L1  *  4  Tl 2  m5  1 


(95) 


and  j  ■  ypl  as  before.  We  define  the  quantity to  be  the  arc¬ 
tangent  of  the  quotient  obtained  by  dividing  the  imaginary  part  of  \ 
by  the  real  part.  Deriving  the  expression  for  c*(t)  by  the  inverse 
Laplace  transform,  we  have  from  Equation  (24) 


ej.(t)  «  0  , 


0  <  t  <  iT 


(96a) 


cx(t) 


e-at. 

09 

Ql  +  %  8^n 

2lTm  x  -  •«-' 

t  ♦  0m  -  mTT 

1 

m»l 

L 

L  L 

J 

e-*aT  sin  [u)t  -  k  <*>  T  ]  ,  (96b) 

£r  ±  t  <  (n  ♦  J)r 


C^t) 


•  I  Bm  |  e'iaT  {8in  [*  “  *  *  T] 

-  e-anT  8in  ♦  /Ba  _  £a)  T  -  «)n  T]}  - 


(96c) 


t  >  (n  ♦  J)T 


Neglecting  the  infinite  series  in  Equation  (96b),  we  obtain  Equations 
(89)  from  Equations  (96)  by  replacing  j  B^  and  m  with  Q2  and  -f 
respectively.  Since  the  model  parameters  were  selected  using  Relation 
(93),  we  must  determine  whether  the  approximations 


are  reasonable. 
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‘9  «  (97) 

q2  **  f^a|  (98) 


Introduce  o-m,  where 
{ 

2a  «)  T2 

4TT2  ra2  ♦  T2(a2  -«)2)J  ’ 

^  4 

For  the  approximation  of  Relation  (9?)  it  may  be  easily  shown  from 
Equation  (95)  that 


(100) 


where  <p  is  given  in  Equation  (88).  let  be  a  non-negative  number 
such  that 


«>  -  ka,  a  .  (101) 

Then  it  is  shown  in  Appendix  A  that 


0  *  °m 


lW  ^  1  (102) 


0  £  <rm  <  IT 


ktf  >  1  .  (103) 


Consequently,  if  the  value  of  u>  is  greater  than  the  bandwidth  of  the 
system,  it  is  necessary  to  establish  a  maximum  for  k^  in  order  that 
the  approximation  of  Relation  (97)  be  reasonable.  By  reasoning  similar 
to  that  preceding  Relation  (93)  we  have  from  Equation  (96c) 
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Ci<t)  «  |Bm  |  e^aT  [l  -  e-*nTl  sin  [a>t  *  /i*  -  i  **]  , 

J  (1C 

t  l  (n  *  i>T  ;  (n  ♦  i)T  »  A  ;  n  T  »  |T  . 

Lei  a  new  quantity  %  be  designated  the  exact  model  phase  lag,  where 
from  Relations  (91 )  and  (104) 

%  -  -♦  -  \yi^  .  dc 


Substituting  the  expression  for  <p  and^XB^  yields 

%  «  *  <rffi  .  (106 

By  the  restrictions  on  the  product  a  T  in  Relation  (50)  we  obtain  from 
Relations  (102)  and  (103) 


0  S  %  <  » 


0  S  %  <  -  T  , 


lCa,  >  1 


Hence  for  W  <;  a  wj  are  assured  that  can  be  no  larger  than  0.271 
radians,  or  approximately  15  degrees.  By  Relation  (108)  this  constraint 
on  the  exact  model  phase  lag  increases  for  >  a.  We  intuitively  assume 
that  a  phase  lag  greater  than  15  degrees  would  not  meet  nominal  engineer¬ 
ing  requirements.  Thus  is  restricted  to  values  less  than  or  equal  to 

the  system  bandwidth.  The  approximation  of  by  -<p  is  then  reason¬ 

able  from  Relations  (100)  and  (102). 

To  Justify  Relation  (98)  Introduce  ,  where 
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(109) 


From  Equation  (95)  it  follows  that 


r  i 

1  MirJ 

(no) 


Let  a  new  tern  y  be  defined  by 


2  2 
<*T2 


(m) 


Then  by  Relation  C 50 ) 


y  < 


i 

161T- 


(112) 


It  is  shown  in  Appendix  B  that 


expi 


k&  -  2k^  -  1 
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Or  U 


&  ■  ll 


£ 

^-14-ijU  6 


(H3) 


ku><  [l+£]  * 


By  Relation  (112) 


[l  *y]*  >  UT 


(114) 


Since  has  already  been  restricted  to  values  less  than  or  equal  to  one^ 
we  have  from  Relation  (113) 
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0.979  <  ~  ^  1.011  ,  0  ^  <  1  (115) 

IT44 

From  Equation  (110)  we  conclude  that  the  approximation  in  Relation  (98) 
is  reasonable. 
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9.  SECOND  ORDER  SYSTEM 

To  note  some  of  the  difficulties  in  selecting  the  delay  line  model 
for  more  complex  systems,  consider  the  transfer  function  of  a  second 
order  system 


0<s) 


b  c 

(s  ♦  b)  (a  +  c) 


b  jta  C 


(116) 


Wiere  b  and  c  are  either  both  real  non-zero  numbers  or  are  non-zero 
complex  conjugates.  Following  a  procedure  similar  to  that  for  Equations 
(9)  through  (19)  we  obtain 


0,(i 


be]; 

e-*We"*T8 

e.  { 

e-  ( n  +i  )bTe-  (n+ J  )T  s 

c  -  bj 

.(b  ♦  s)  J 

u 

fi 

e-^e-ift8  -  e-(n4)cTe-(n^)TB  ' 
",  t  if  (  *2(c  ♦  o)2  i 

<c * *>  1  [‘‘TFirJ  J 

Let  the  quantities  g£[t]  and  gj[t]  be  given  by 

«:w  ■  rttst  [!-•■“] 


(117) 


(118) 
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s*  w  •  t1  -  •'“]  •  <u,) 

By  omitting  the  infinite  products  in  Equation  (11?)  it  may  be  shown  that 
the  approximate  model  response  to  a  unit  step  input  becomes 

ciA^)  "  «"^bT  «b  [*  ~  £T]  «(*  "  £T) 

_  a-(n4)bT  g*  [t  -  (n  ♦  £)t]  u[t  -  (n  ♦  $T )] 

(120) 

-  e"^cT  &Z  [t  -  *T]  u(t  -  |T) 

♦  [t  -  (n  ♦  **]  u[t  -  (n  ♦  J)t] 


The  final  value,  V,  of  the  approximate  model  response  is  given  by 


The  problem  of  selecting  T  to  maximize  V  is  more  complicated  than  for  the 
first  order  system.  Taking  the  derivative  of  the  expression  for  V  with 
respect  to  T,  equating  the  result  to  zero,  and  rearranging,  we  obtain 

.  LzJ* LLA)  »-c°t  .  (122) 

e-^cT  1  -  (2n  ♦  1)  e-taT 

An  explicit  expression  for  T0  in  terms  of  n,  c,  and  b  cannot  in  general 
be  found  from  Equation  (122).  If  the  input  to  the  system  is  a  ramp,  the 
final  slope  S  of  the  approximate  model  response  is  also  given  by  the  ex¬ 
pression  for  7.  Further  study  is  therefore  necessary  to  find  a  method 
for  selecting  the  model  parameters. 


56 


10.  SUMMARY  AND  CONCLUSIONS 

A  delay  line  model  of  a  linear  system  is  selected  such  that  the 
responses  of  the  model  to  comnonly  occurring  inputs  are  closest  to 
corresponding  responses  of  the  system.  It  is  found  for  each  input  that 
the  number  n  of  model  delay  elements  should  be  chosen  as  large  as  prac¬ 
tical  for  applications.  The  selected  value,  Tc,  of  the  delay  time,  T, 
is  then  given  by 

Ln  (2n0  ♦  1) 

T  .  - 

0  a  no 

where  a  is  the  system  bandwidth  and  no  the  selected  n.  For  the  periodic 
input  studied  the  signal  frequency  is  limited  to  values  equal  to  or  less 
than  the  system  bandwidth  in  order  that  the  model  be  acceptable. 

The  expression  for  T0  is  obtained  by  a  technique  using  infinite 
products  previously  employed  successfully  by  Oldenburger  and  Goodson^ 
in  distributed  parameter  studies.  ?y  assuming  a  small  value  of  the 
model  delay  time  the  infinite  product  appearing  in  the  model  transfer 
function  is  neglected  and  a  simplified  model  representation  obtained. 

It  is  shown  that  this  approximation  is  reasonable  for  the  transient 
inputs  if  aT  <  The  same  is  true  for  the  periodic  input  if  in  ad¬ 
dition  tne  input  frequency  is  less  than  or  equal  to  the  system  band¬ 
width.  Although  the  appro xima t ion  method  is  quite  useful  in  the  first- 
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order  system  study,  new  difficulties  noted  at  the  conclusion  of  the 
investigation  arise  for  more  complex  systems.  Since  T0  is  a  function  of 
the  system  parameters  as  well  as  the  selected  number  of  delay  elements, 
the  difficulty  in  model  selection  increases  as  the  number  of  system 
parameters  increases. 
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APFENDIX  A 

ESTABLISHMENT  OF  CONSTRAINTS  ON  crm 


Let  the  quantity  be  given  by 

2a  u>T2 _ 

%  "  AT  2  m2  ♦  t2(*2  -  <£) 

Then  from  Equation  (99) 

ao 

o-»  ■  ZZ  tan-1  (xJ 

in»l 

From  Equation  (95)  it  is  easily  shown  that 

0  *  <rB  <  T  ,  ^  ° 

for  finite  Hereafter  we  consider  only  finite,  non-negative 

of  a).  By  Equation  (123)  it  follows  that 


-1  £  %  £  0  ,  1  i  m  i  - 

2J 


csfi  -  2a«)  -  a2 


T  f  o  -  2l 

*  t  r 

—  £  *m  <  -1  .  -gjdtiT  ~2a a)  -  a^ 

<  m  <  — 

*  2TL 

1  <  Xm  < 


■  ^2- 


1  T 

*  <£  Bl  <  —— 

zr 


a)2  ♦  2a  a)  - 


(123) 


(124) 


(125) 

values 


(126) 


(127) 

a2]* 

J (128) 
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0  £  *ra  £  1 


T  2 

s®2 


2a  a) 


(129) 


We  recall  that  m  in  the  above  equations  is  a  positive  integer.  Let  mi 
be  the  first  positive  integer  (FPl)  less  than  or  equal  to,  depending  on 
the  value  of  a),  the  upper  limit  on  m  in  Relation  (126).  From  Relation 
(127)  let  mg  be  the  FPI  greater  than  this  expression.  The  term  m3  is 
the  FPI  less  than  or  equal  to  the  upper  limit  on  m  in  Relation  (128), 
and  so  on  with  the  next  two  relations  through  m^.  With  the  series 
representations 


tan“^  x 


.  .  • ,  ^  |  I 


•  X  £  tan_1  *  < 

4  a  **  4 


(130) 

51T 


tan-1  x 


*T  4- 


(131) 


X-  ^  tan"1  x  <  HL 
4  *  =4 


tan"1  x 


1  3  .  1  5 

3  x  5  x 


x2"  js.  1 


-  X  £  tan'1  £  X 


(132) 


we  keep  only  part  of  each  series  to  give  from  Equations  (126)  through 
(129) 
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0  £  <r 


£  < 


la  [  IT 4  ** "  3  *2]  *  |^2 


^*4  [iT  ‘  i  *  3*jf]  *  ^  W  j  * 


(133) 


If  u)  £  a,  then  m  is  defined  only  for  m  £  m£.  For  instance,  in  Relation 
(129)  let  u)  ■  k^j  a,  where  0  £  ku)  g  1,  Then  since  a  T  <  £ 


zH*2  * 2au;  -•*!*<  “ij 


Hfence  from  Relation  (133) 


0  &  <rm  < 


M  »  ^  ^  a  . 


(134) 


Let  quantities  P  and  Q  be  given  by 


P  - 


2a  a)  T8 

2TZ 

T2(a2  -  <Q2) 

4T2 


(135) 

(136) 


Then  from  Equation  (123) 


Xm  "  P 


Li. 

[m2  ♦ 


(137) 


It  follows  that 
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XlB  <.  P 


^  a 


By  Equation  (55)  we  obtain 


[?] 


■■■6 

Therefore  from  Relation  (50) 


1  u> 
12 


tO  i  a 


0  s  <r«  “  Is 


Cl) 


£  * 


and  Relation  (102)  is  established. 

For  oi  >  a  we  have  by  Relation  (125) 


<rm 


T  . 


ku)  >  1 


(138) 


(139) 


(340) 


(141) 


This  establishes  Relation  (103) 
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APPENDIX  B 

ESTABLISHMENT  OF  CONSTRAINTS  ON 


Employing  Equations  (101)  and  (ill)  we  obtain  from  Equation  (109) 


(142) 


It  follows  that 


X  4  X  *  l*>^l 


(143) 


Let  k-^  be  defined  by 


*1 


■m 


(144) 


Then  for  the  quantity  in  the  first  pair  of  vertical  bars  in  Relation 

(143), 


>  0  ,  0  S  k<o 


<  k! 


(145) 


Fo£  tjne  remainder  o£  this  work  the  limits  on  in  Equation  (145)  are 
assumed. 
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Introduce  Qjj,  where 

y 

<3m  “  -2"  [kS-  2k#)-  l]  .  (146) 

Whence,  by  Relation  (145 )> 

H  *  f  [l  -  Si]  •  (W) 

From  Equation  (146)  and  Relation  (112), 

“  gfST  <  ^  <  1  (U8) 

for  any  m.  Let  the  quantity  Ug,  be  given  by 

%  "  £  141  l1  -  Sn]  •  (149) 

m*! 

Then 

••  . 

If  [l  -  9m  ]  -  «Um  .  (150) 

m-1  L 

The  following  established  relationships  are  noted*32 

In  (1  ♦  x)  g  x  ,  x  £  0  (151) 

— ~  £  In  (1  -  y)  £  -y  ,  0  <  y  <  1  .  (152) 

y  *  ^ 
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Let  vm  be  defined  by 


Hence 


T  [I  -  R»]  ■  *Vm 

a«l 

Define  the  quantity  R  by 

R  ■  y|kS  -  l| 


(160) 


(161) 


Thus 


-ISnl  .  * 

|Rm|  -  1  R  -  m2 


(162) 


Since  from  Relation  (112) 

y  1 4  -  i|  *•  m2 


(163) 


for  all  m,  we  have 


12 


(164) 


We  note  that 
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We  substitute  the  expression  for  R  in  Relation  (166)  and  employ  Elation 
(55 )«  Then  by  Equation  (l60) 


Relations  (155)  and  (167 )  establish  Relation  (113). 


